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Abstract 

We present a study of the tensorial structure of the hadronic matrix elements of the angular 
momentum operators J. Well known results in the literature are shown to be incorrect, and we have 
taken pains to derive the correct expressions in three different ways, two involving explicit physical 
wave packets and the third, totally independent, based upon the rotational properties of the state 
vectors. Surprisingly it turns out that the results are very sensitive to the type of relativistic spin 
state used to describe the motion of the particle i.e. whether a canonical (i.e. boost) state or a 
helicity state is utilized. We present results for the matrix elements of the angular momentum 
operators, valid in an arbitrary Lorentz frame, both for helicity states and canonical states. 

These results are relevant for the construction of angular momentum sum rules, relating the 
angular momentum of a nucleon to the spin and orbital angular momentum of its constituents. It 
turns out that it is necessary to distinguish carefully whether the motion of the partons is character- 
ized via canonical or helicity spin states. Fortunately, for the simple parton model interpretation, 
when the proton moves along OZ, our results for the sum rule based upon the matrix elements of 
J z agree with the often used sum rule found in the literature. But for the components J x , J y the 
results are different and lead to a new and very intuitive sum rule for transverse polarization. 



1. INTRODUCTION AND SUMMARY OF RESULTS 

Sum rules, relating the total angular momentum of a nucleon to the spin and orbital 
angular momentum carried by its constituents, are interesting and important in under- 
standing the internal structure of the nucleon. Indeed it is arguable that the main stimulus 
for the tremendous present day experimental activity in the field of spin-dependent struc- 
ture functions was the surprising result of the European Muon Collaborations polarized DIS 
experiment in 1988 [1], which, via such sum rules, led to what was called a "spin crisis in 
the parton model" [2], namely the discovery that the spins of its quarks provide a very small 
contribution to the angular momentum of the proton. A key element in deriving such sum 
rules is a precise knowledge of the tensorial structure of the expectation values of the angular 
momentum operators J, in a state \p , a) of the nucleon, labeled by its momentum p, and 
with some kind of specification of its spin state, denoted here non-commitally by o . 

In a much cited paper [3], Jaffe and Manohar stressed the subtleties involved in deriving 
general angular momentum sum rules. As they point out, too naive an approach leads 
immediately to highly ambiguous divergent integrals, and a careful limiting procedure has 
to be introduced in order to obtain physically meaningful results. In this it is essential to 
work with non-diagonal matrix elements (p',cr\J \p,cr) and, as discussed below, this can 
have some unexpected consequences. Jaffe and Manohar comment that to justify rigorously 
the steps in such a procedure requires the use of normalizable wave packets, though they do 
not do this explicitly in their paper. 

In a later paper [4], Shore and White utilized the approach of Ref. [3], including an 
explicit treatment with wave-packets, to derive some far reaching conclusions about the role 
of the axial anomaly in these sum rules. 

We shall argue that despite all the care and attention to subtleties, there are flaws in 
the analysis in [3] and the results presented there are not entirely general. Indeed there are 
cases where the results of [3] are incorrect. This, in turn, throws doubt upon some of the 
conclusions reached in [4], which we will examine. 

The bulk of our analysis is based on a straightforward wave-packet approach. However, as 
we explain, this is rather subtle for particles with non-zero spin. The key points underlying 
our results are: 

(1) Our wave packets are strictly physical, i.e. a superposition of physical plane-wave 



states. This requirement turns out to be incompatible with some of the steps in [3]. 

(2) We give a careful treatment of the Lorentz covariance properties of the matrix 
elements involved in the subsidiary steps of the anaylsis. This leads to tensorial structures 
which differ in some cases from those in [3]. 

(3) Because our results differ from [3] we have looked for and found a totally independent 
way to check our results. This does not use wave packets and is based upon the transfor- 
mation properties of momentum states under rotations. This very direct approach holds for 
arbitrary spin, whereas in the wave packet treatment we are only able to deal with spin | 
particles. It also brings to light some peculiar and unintuitive properties of helicity states, 
and this must be taken into account when deriving spin sum rules. This is important since 
we have to deal with gluons in our sum rules. 

Our results for the matrix elements of J are as follows: 

For a massive particle of spin | with 4-momentum p in a canonical spin state \p, s) (i.e. 
in a 'boost' state of the kind generally used in textbooks on Field Theory e.g. in Bjorken 
and Drell [5], or in Peskin and Schroeder [6] ) , where s/2 is the spin eigenvector in the rest 
frame (s 2 = 1), we show that, for the forward matrix elements, 



- St + i(p X V p )i 



8\p' - p) . (1.1) 



(p',s\J l \p,s) = 2p°(2n) 3 

The states \p, s) are normalized conventionally to 

(p',s\p,s)=2p°(27f) 3 5 3 (p'-p) (1.2) 

and we note that in the rest frame 

\0,s) = V 1 r £ /2 (R(s))\0,m) (1.3) 

where | 0, m) has spin projection m along the ^-direction in the rest frame and R(s) rotates 
a unit vector in the z-direction into s by first a rotation about y and then a rotation about 
z. 

For the purpose of deriving sum rules our result for the matrix elements non-diagonal in 
the spin label is actually more useful, namely, for a spin \ particle 



{p',m'\Ji\p,m) = 2p (2wy 



1 d 



5^\p'-p). (1.4) 



The generalization of these results for arbitrary spin is given in Eq. (6.13). 

Helicity states are more suitable for massless particles such as gluons. Using the Jacob- 
Wick conventions for helicity states [7] we find a surprisingly different result, namely, 

(p', X'\Ji\p, A) = (27T) 3 2p [Xin(p) + i(p x V P )J S 3 (p' - p)5 xx , (1.5) 

where 

i] x = cos(0) tan(6*/2), r\ y = sin(0) tan(6 l /2), r\ z = 1. (1.6) 

and (9, <fi) are the polar angles of p . 

The first term in Eq. (1.1) differs from the results of Jaffe and Manohar [3] . If we rewrite 
their expression (see Eq. (4.23)) in terms of the independent vectors p and s, we find, for 
the expectation value 



«>■>« = !£■?(<%* - «^ - -f^W (P- *•) M 

to be compared to 

(Ji) = \si (1.8) 

arising from the first term in Eq. (1.1). In general these are different. However, one may 
easily check that if s = p the Jaffe-Manohar value agrees with Eq.(1.8), while if s J_ p they 
are not the same. 

The agreement for s = p is consistent with the much used and intuitive sum rule 
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AS + AG + (L q ) + (L G ) (1.9) 



In the case that slpwe find a new sum rule. For a proton with transverse spin vector St 
we find 

3 = 3 E JdxA T q a (x) + £ (L ST > a (1.10) 

where L ST is the component of L along s-p. The structure functions Aj-g a (x) = /i?(a;) are 
known as the quark transversity or transverse spin distributions in the nucleon. Note that no 
such parton model sum rule is possible with the Jaffe-Manohar formula because, as p — > oo, 
Eq. (1.7) for i — x,y diverges. 

The result Eq. (1.10) has a very intuitive appearance, very similar to Eq. (1.9). 

The organization of our paper is as follows: In Section 2 we explain why the calculation 
of the angular momentum matrix elements is so problematical. Because of the unexpected 



sensitivity of the matrix elements of J to the type of spin state used, and because we are 
forced to use helicity states for gluons, Section 3 presents a resume of the difference between 
Jacob- Wick helicity states and canonical ( i.e as in Bjorken-Drell) spin states. Further, 
given that our results disagree with one of the classic papers in the literature, we have felt it 
incumbent to summarize, in Section 4, the treatments of Jaffe-Manohar and Shore- White, 
pointing out the incorrect steps in these derivations. 

In Section 5 we present a detailed wave-packet derivation of the structure of the matrix 
elements of J for spin |, first for a relativistic Dirac particle, then in a field theoretic 
treatment. We comment here on the claims made in Shore- White on the role of the axial 
anomaly in the structure of the matrix elements. Sections 4 and 5 are heavy going, and the 
reader only interested in a quick and direct derivation of the key results should skip these 
and read Sections 6, 6.1, 6.3 and 7. 

In Section 6 we confirm the results of Section 5 in a completely independent approach, 
which is valid for arbitrary spin, based on the rotational properties of canonical and helicity 
spin states. We also prove that our results are in conformity with the demands of Lorentz 
invariance. 

In Section 7 we derive the most general form of an angular momentum sum rule for a 
nucleon and show that it reduces to the standard, intuitive, sum rule for J z when the nucleon 
is moving along OZ . We also derive a new sum rule for a transversely polarized nucleon. 

2. THE ORIGIN OF THE PROBLEM 

In the standard approach one relates the matrix elements of the angular momentum 
operators to those of the energy-momentum tensor. 

Let T^ v (x) be the total energy- momentum density which is conserved 

d,T^(x) = 0. (2.1) 

Later we shall distinguish between the conserved canonical energy- momentum tensor Tq V , 
which emerges from Noether's theorem, and which is, generally, not symmetric under \x <-> is, 
and the symmetrised Belinfante tensor T MJ/ , which for QCD is given by 

^(x) = i(^(x) + T c ») (2.2) 

and which is also conserved. For the moment, however, this distinction is irrelevant. 
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Being a local operator, T^ u (x) transforms under translations as follows 

T^(x) = e ix - p T^ u (0)e- lx - p , (2.3) 

where the P^ are the total momentum operators of the theory. 

By contrast the various angular momentum density operators which are of interest, the 
orbital angular momentum densities 

<?(*) ee x»T»\x) - x x T^(x) (2.4) 

or the version constructed using the symmetrised stress-energy tensor, 

M" u \x) = x v T^(x) - x x T^(x) (2.5) 

are not local operators (we shall call them compound) and do not transform according to 
Eq. (2.3). 

Note that, strictly speaking, the operators relevant to the angular momentum are the 
components M ^ where i,j are spatial indices. However, for reasons of simplicity in utilising 
the Lorentz invariance of theory, the authors of Ref. [3] prefer to deal covariantly with the 
entire tensor M^ vX . We shall loosely refer to them also as angular momentum densities. 

The total angular momentum density is 

J^\x) = M^\x) + MjU(x), (2.6) 

where the structure of M^ in depends on the type of fields involved. From Noether's theorem 
J fluX (x) is a set of conserved densities, i.e., 

d^ v \x) = 0. (2.7) 

As a consequence of the densities being conserved, it follows that the total momentum 
operators 

F- S j **<*(*) (2.8) 

and the total angular momentum operators J, 

J z = J 3 = J 12 , (cyclical) (2.9) 

with 

J ij = [ d 3 xJ 0t:i (x) (2.10) 



are conserved quantities, independent of time. 

The relationship between the M> IuX (x), constructed using the symmetrical energy- 
momentum density and the J^ vX (x) constructed from the canonical energy-momentum ten- 
sor is extremely interesting and will be commented on later in Section 4.2. One can show 
(see e.g. Ref. [4]) that 

M 0ij (x) = J 0ij (x) + [E of M terms] + [divergence terms]. (2.11) 

The [E of M terms] vanish if it is permissible to use the equations of motion of the theory. 
The [divergence terms] are of the form d a F a0lJ (x) . 

As mentioned we shall be primarily interested in the expectation values of the physical 
operators, i.e. in their forward matrix elements. If F a0 ^(x) were a local operator, it would 
follow directly that the forward, momentum-space, matrix elements of the divergence terms 
in Eq. (2.11) vanish. But it is not a local operator. Nonetheless, a careful treatment using 
wave packets [4] demonstrates that the forward matrix elements do indeed vanish. See 
Section 4 below. 

Dropping, as is customary, the [E of M terms], we shall thus assume the validity of 



(p,a\ / d 3 xM 0ij (x,0)\p,a} = (p,a\ / d 3 xJ 0ij (x, 0)\p, a). (2.12) 

We shall return to this question in Sections 5.1 and 5.2. 

Of primary interest are the matrix elements of the angular momentum operators J k or, 
equivalently, the J l K Consider the forward matrix element, at t — 0, 

M 0ij (p,s) = (p,s\ ( d 3 xM 0ij (x,0)\p,s) (2.13) 

= f d 3 xx l (p, s\e ip - x T^ (0)e- lP - x \p, s) - (i <-► j) 

= I d 3 xx i (p,s\T 0j \p,s)-(i^j). (2.14) 

The integral in Eq. (2.14) is totally ambiguous, being either infinite or, by symmetry, zero. 

The essential problem is to obtain a sensible physical expression, in terms of p and s, 

for the above matrix element. The fundamental idea is to work with a non-forward matrix 

element and then to try to approach the forward limit. This is similar to what is usually 



done when dealing with non-normalizable plane wave states and it requires the use of wave 
packets for a rigorous justification. 

It will turn out that the results are sensitive to the type of relativistic spin state employed, 
so in Section 3 we present a brief resume of the properties of relativistic spin states. We 
then proceed to discuss the approaches of [3] and [4] in Section 4, where we shall comment 
on the dubious steps in these treatments. The most crucial error in these treatments is the 
mishandling of the matrix elements of a covariant tensor operator. If T^ A transforms as 
a second-rank tensor its non-forward matrix elements do not transform covariantly. This 
was the motivation, decades ago, for Stapp to introduce M- functions [11]. Namely, the 
covariance is spoilt, for canonical spin states by the Wigner rotation, and, for helicity states 
by the analogous Wick helicity rotation [10]. Only by first factoring out the wave- functions 
(in our case Dirac spinors) i.e. by writing 

(p',S / \T^\p,S)=u(p',S')T^(p',p)u(p,S). (2.15) 

does the remaining M-function, in this case T ftv {p' ) p) ) transform covariantly. For local 
operators the transformations of the spinors u and u cancel between themselves for forward 
matrix elements and so the result does have the naively expected tensor expansion. This is 
not true in general for compound operators, in particular the angular momentum and boost 
operators. 

3. RELATIVISTIC SPIN STATES 

The definition of a spin state for a particle in motion, in a relativistic theory, is non- 
trivial, and is convention dependent. Namely, starting with the states of a particle at rest, 
which we shall denote by [0, m), where m is the spin projection in the z-direction, one defines 
states \p, a) for a particle with four-momentum p by acting on the rest frame states with 
various boosts and rotations, and the choice of these is convention-dependent. The states 
are on-shell so p 2 = M 2 

There are three conventions in general use [8] 

a) Canonical or boost states as used e.g. in Bjorken and Drell [5] or Peskin and Schroeder 

[6] 

\p,m) = B(v)\0,m) (3.1) 



where B(v) is a pure boost along v = p/po, and p = (p, 8, <fi) denotes the three-vector part 
of p^ 

b) Jacob- Wick helicity states [7] 

\p, \) JW = RM R y (9) R z {-4>) B z (v)\0, m = \) (3.2) 

where B z is a boost along OZ, and the later introduced, somewhat simpler 

c) Wick helicity states [10] 

\p, A) = R z (<j>) R y (9) B z (v)\0, m = \) (3.3) 

From the canonical states of spin-| one can construct the states 

\p,s) = B(v)\0,s) = B(v)V^ /2 (R(s))\0 : m) (3.4) 

which, in the rest frame, are eigenstates of spin with spin eigenvector along the unit vector 
s. (The rotation R(s) was explained after Eq. (1.3)). 

The canonical states, with their reference to a rest frame, are clearly not suitable for 
massless particles like gluons. Helicity states, on the other hand, can be used for both 
massive and massless particles. However, it turns out that the results for the canonical 
states are much more intuitive, so we will generally use them for M^O. 

The reason we are emphasizing this distinction between canonical and helicity states is 
that the matrix elements of the angular momentum operators between helicity states are 
quite bizarre! Since, for arbitrary p, helicity states are just linear superpositions of canonical 
states, one may wonder why this is so. It results from the facts i) that the coefficients in the 
linear superposition are p-dependent, i.e. depend upon the polar angles of p and ii) that 
the matrix elements of the angular momentum operators contain derivatives of (^-functions, 
and these, as usual, must be interpreted in the sense of partial integration, i.e. 

f(p,p')^S 3 (p-p f ) = S 3 (p-P')^f(p,p') (3.5) 

dpi dpi 

The technical details are explained in Section 6. 

In almost all studies of hard processes, where a mixture of perturbative and non- 
perturbative QCD occurs, nucleons are taken to be in helicity states moving with high 
energy along the z-axis, and typically one is utilizing matrix elements of local products of 
quark or gluon field operators between these states. For these operators there is no problem 
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in dealing with diagonal matrix elements. But when it comes to an angular momentum sum 
rule for the nucleon , care must be taken to decide whether one is dealing with helicity states 
\p z ,X), where p z = (E,0,0,p) or with canonical states \p z ,s z ), where s z = (0,0, 2A). The 
point is that even though the initial states are the same, 

\p z , X) = \Pz, s z ) (3.6) 

the singular nature of Jj forces one to deal with non-diagonal matrix elements i.e. to utilize 
(p', a\ where p' is not along the z-axis, and for these 

(p',X\^(p',s z \ (3.7) 

In this paper we show that it is possible to give a rigorous derivation of the structure of 
the expectation values for canonical states 

(p, s\Ji\p, s) = C p '^ p (p', s\Ji\p, s) (3.8) 

where s is a unit vector along the rest frame spin eigenvector, and for helicity states 

(p, X\Ji\p, X) = C p ^ p (p', X\Ji\p, A). (3.9) 

In general, for the arbitrary i th component of J, for spin | 

(p, \\Ji\p, X) ^(p,s = 2Xp\Ji\p, s = 2Xp) (3.10) 

even though s lies along the direction of p in both cases, and even if p is along OZ where 
Eq.(3.6) holds. Only for the specific component of J along p do the matrix elements agree, 
i.e. for arbitrary p, 

(p,X\J- P\P, X) = ( P ,s = 2Xp\J • p\p, s = 2Xp). (3.11) 

In using the sum rules based on Eq.(3.8) or Eq.(3.9) for arbitrary % to test any model of the 
nucleon in terms of its constituents, it is essential to construct wave-functions appropriate 
to the type of spin state being used for the nucleon. The equations Eq.(3.9) and Eq.(3.8) 
contain delta functions and derivatives of delta-functions and this is the reason for the special 
care required. Throughout this paper, with the exception of the discussion in Section 6 and 
in Section 7, we will use canonical spin states. In these latter sections we shall utilize 
Jacob- Wick helicity states for the massless gluons. 

In Section 4 we summarize the treatments of Jaffe-Manohar and Shore- White. The 
approaches in Refs. [3] and [4] are different, so we shall present both in some detail and will 
comment upon the dubious steps. 

11 



4. THE JAFFE-MANOHAR AND SHORE- WHITE APPROACHES 

These authors employ the standard approach of trying to relate the matrix elements of 
the angular momentum operators to those of the energy-momentum density operator, and 
utilize the version Eq. (2.5) based on the symmetrized energy-momentum tensor. 

4-1. The Jaffe-Manohar treatment 

In order to make efficient use of the Lorentz invariance the authors of [3] prefer to label 
their states using the covariant spin 4-vector S and to consider the entire tensor J^ vX and 
to integrate the tensor densities over four-dimensional space time, i.e. they consider 

M^ x {p,k,p',S)= f d 4 xe ikfX {p',S' = S\M> M/X {x)\p,S) (4.1) 

and eventually take the limit /c M — > 0. In Ref. [3] the LHS is written in the abbreviated form 
Ai^ uX (p, k, s), but we shall use the above notation for clarity. Note that Jaffe and Manohar 
use the notation s to mean the covariant spin which we denote by S. It is IjM times the 
expectation value of the Pauli-Lubanski operator; see e.g. [8]: 

W„ = -\^ vpa P v J pa (4.2) 

and 

(2n) 3 p M5 3 (p' - p)S, = (p', s\W,\p, s) (4.3) 

In terms of components 

p • s 



S, 







M 



S l = st + p* P ' S (4.4) 

F M(p + M) v ' 

Also we take S 2 = — 1 while Ref. [4] takes the covariant normalization to be — M 2 . 
The Lorentz invariant normalization of the states is conventional 

(p',S\p,S) = (27r) 3 2p 5 3 (p f -p). (4.5) 

The extra integration in Eq. (4.1) J dt, is argued to be harmless, leading to an infinite 5(0) 
which cancels out when calculating a genuine expectation value. 
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Comment 1 It will be seen in Section 5.3 that the choice S' = S as done in [3] is not 
consistent in a proper wave-packet treatment. 

Analogously to the steps leading to Eq. (2.14) we have 

M" uX (p, k,p',S) = f d 4 xe ix < k ~ p+p,) x u (p', S\T^(0)\p,S) - (u <-► A) 

. d 



d x 



-i(2ir) 



D ix-(k—p+p') 



' dk v 



(p',S\T» x (0)\p,S)-(is~\) 



dk„ 



[5\k-p + p f )(p',S\T^(0)\p,S)]-(u^\). (4.6) 



Comment 2 The last step in Eq. (4.6) can only be justified if p' in Eq. (4.1) is considered an 
independent variable. We may not take p' = p — k. Once this is recognized it is no longer 
so evident that Eq. (4.1) followed by the limit k^ — ► provides a natural definition of the 
ambiguous forward matrix element. 
Continuing from Eq. (4.6) one has 



M^ x (p,k,p',S) 



-z(27r) 



(p-k,S\T»\0)\p,S)^5 4 (k-p + p') 



d 



+ S\ k -p + p') — {p-k,S\T^\p,S) 



(!/<-> A). (4.7) 



In Ref. [3] the limit k^ — » is given as 



M^\p,0,p',S) 



-i(2v) 



(p,S\T»\0)\p,S)d»5\0) 



d 



+S\0) w {p-k,S\T^\p,S) 



-(i/<-» A). 



(4i 



This form is a little puzzling, given that p' ^ p — k in Eq. (4.1), as discussed in Comment 2. 
We thus prefer to write Eq. (4.8) in the form 



M^\p, 0,p',S) = -i(2ir) 



(p-k,S\T» x (0)\p,S) lim ^-5\k-p + p') 

fcM^o ak v 



d 



+S\p'-p)—(p-k,S\T^\p,S) 



^-►oJ 



(i/<-*A). (4.9) 



The highly singular first term in Eq. (4.9) can only be understood in a wave-packet 
analysis. It corresponds to the angular momentum about the origin arising from the motion 
of the center of mass of the wave-packet, and has nothing to do with the internal structure 
of the nucleon. Thus we will take as the definition of the ambiguous forward matrix element 
in the Jaffe-Manohar (JM) approach 



M^\p,0,p',S) = -i(2ir) 4 5(0)5 3 (p' -p) ^{p-k,S\T» x (0)\p,S) 
-(*<-> A), 



fc'^O 



(4.10) 
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where we have used the fact that for the on mass-shell momenta p' = p forces p' = po- 
The last part of the analysis concerns the structure of the matrix element of T /iA (0). For 

Eq. (4.10) we require an expansion in k^ of (p — /c,5|T MA (0)|p, S) up to terms linear in k^ '. 

It is at this point that the choice in Eq. (4.1) of S' = S becomes significant: it greatly 

simplifies the tensorial structure of the expansion. 
Following Ref. [3] one writes, with P M = p M + \k^, 

(p - k, S\T» x (0)\p, S) = A {k 2 )P^P x + iA 1 {k 2 )[e^ a P x + e Aa/3a P»}k a Pp S a + 0{k 2 ), (4.11) 

where A and A\ are scalar functions of k 2 . The terms in Eq. (4.11) are chosen to respect 
the relation 

k,(p-k,S\T^(0)\p,S) = 0. (4.12) 

(To see that k ■ P = one should recall that the nucleon states are physical states with 
(p-A;) 2 =j9 2 = M 2 .) 

For the forward matrix element 

(p,S\T^(0)\p,S)=A (0)p»p x . (4.13) 

Comment 3 While Eq.(4.12) is correct the crucial expansion (4.11) is not. The reason 
is the following. The RHS of Eq. (4.11) has been constructed to transform under Lorentz 
transformations as a genuine second-rank tensor on the grounds that T mA (0) transforms as 
a second-rank tensor. But the non-forward matrix elements of a tensor operator do not 
transform covariantly, as was explained in the discussion of Eq. (2.15). We shall see the 
consequences of this in Section 5.1. 

Continuing with the derivation in Ref. [3], we have from Eqs. (2.8) and (4.5) 

(p',S\Jd 3 xT 00 (x)\p,S) = (p',S\H\p,S) 

= 2p 2 (2n) 3 5 3 (p'-p), (4.14) 

where we have used the fact that P° is the total energy or Hamiltonian operator. But the 
LHS, taking t = 0, equals 

f d 3 xe ix < p - p '\p',S\T 00 (0)\p,S} = (2tt) 3 5 3 (p'-p)A (0)p 2 (4.15) 

from Eq. (4.13). 

14 



Comparing Eqs. (4.14) and (4.15) yields 

A (0) = 2. (4.16) 

Now one uses Eq. (4.11) to calculate the derivative needed in Eq. (4.10). The result in 
Ref. [3] is 

M^ x (p,0,p',S) = (2tt) 4 5(0)5 3 (j/ -p)Ai(O) \2p^t Xv ^ - ft^ u + pV^ ff ] p^S,. (4.17) 

Lastly the value of Ai(0) is found by choosing a nucleon state at rest and spin along OZ 
This is an eigenstate of J z 

J 2 |0,£> = ±|0,£> (4.18) 

Then with S= (0,0,0,1) 



(p',z\ I d 4 xM 012 (x)\0,z) = J dt(p',S\J z \0,S) 



= \ {2n) 4 2M5 3 {p')5(0). (4.19) 

But from Eq. (4.17) the LHS is just equal to "I" 

M 012 {0,0,p',S) = (2tt) 4 5(0HV)^i(0)2M 2 (4.20) 

so 

AU0) = . (4.21) 

K ' 2M K J 

Finally, then, in the Jaffe-Manohar treatment the interpretation given to the forward 
matrix element of the angular momentum operator is 



(p',S\ I d 4 xM^\x)\p,S) = M^ x (p,0,p',S) 

(2ir) 4 5(0)5 3 (p'-p) 



[2p^e Xu(3a - p v e ^ Xf3a + p x e^ a ] ppS a . 

(4.22) 

Equation (4.22) is meant to provide a general basis for angular momentum sum rules. 
So, for example, if we have a theory of the nucleon in terms of quark and gluon fields and we 
construct the operator M fluX from these fields, then the requirement that our M^ vX satisfy 



T There is a factor of M 2 missing in the value of A\ (0) given in [3] . 
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Eq. (4.22) for an arbitrary state of the nucleoli yields a set of conditions on some of the 
elements of the theory. 

As indicated in the comments, there are flaws in the derivation and Eq. (4.22) is incorrect. 
We shall present the correct result, in a wave-packet treatment in Sections 5.1 and 5.2. There 
our states or wave-functions will be normalized to 1 so that we calculate actual expectation 
values. Moreover, as will be explained in Section 4.2, the wave-packet approach seems 
only able to deal with the physically relevant operators J d 3 xM 0lJ (x) . For these conserved 
densities the J dt in Eq. (4.1) is simply equivalent to the factor 2ir5(0) in Eq. (4.17). Thus, 
dividing the expression (4.22) by this factor and by the normalization given in Eq. (4.5) we 
have, for the expectation values in the JM-treatment 

(p,S\fd 3 xM 0i i(x)\p,S) 



= -JIT- W e3l<i ° ~ P te ° 3P<T + P 3 e 0tf3cT ] p p S a (4.23) 

jm 4M Po 



(p, S\p, S) 

which can be compared directly with the result we shall obtain in Section 5.1 and Section 5.2. 
In Section 6 we shall give a completely independent corroboration of these results. 

We turn now to the treatment of Shore and White, which basically follows the pattern of 
Ref . [3] , but attempts to put the argument on a rigorous footing via the use of wave-packets. 

4-2. The Shore- White treatment 

As already mentioned, the authors of Ref. [3] remark that a wave-packet approach is 
needed to justify some of the manipulations involved, more precisely, to get rid of the 
unwelcome derivatives of 5-functions. This is done in Ref. [4], but, as we shall see, the 
treatment also suffers from some of the incorrect elements commented on in Section 4.1. 
The notation in Ref. [4] differs somewhat from that of Ref. [3], so we will rephrase the 
notation in Ref. [4] to match as closely as possible the development in Ref. [3]. 

The authors of Ref. [4] try to give a sensible definition to the forward matrix elements 
defined in Eq. (2.14), i.e. 

M 0ij (p, S) = (p, S\ f d 3 x [x'T^ix) - x j T 0i (x)] \p, S) (4.24) 

by utilizing a wave packet. Thus they define 

J (27T)V 2( ?0 
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where drops rapidly to zero as \q — p\ — ► oo, and they interpret Eq. (4.24), modulo 
normalization, as 



M%i(p,S) = (<f>(p,S)\ / d 3 x [x^W-^(x)] \p,S). 



(4.26) 

It turns out to be sufficient to use just one wave packet, either for the inital or the final 
state. 

Note that this differs from J\4 /J,I/X (p,0,p',S) in two respects. Firstly, consideration is 
given here only to the spatial elements i, j of the tensor and secondly, the integral is over 
three-dimensional space. The latter difference is not significant given that the operators in 
Eq. (4.24) are supposed to be time-independent. 

The derivation then runs as follows: 

73, 



M%(p,S) = j ^<P[( q - P Y][x^<^( q ,S\T^(0)\p,S) 



[i^j 



Iwf^-^I^ dr^' (p ~ 9) ) <^i T ° j '(°)b' 5 > 



-t* <-► j)- 



(4.27) 



Exchanging the order of integration, then integrating by parts w.r.t. q, and discarding the 
surface terms at g» = ±oo, 

M°^(p,S) = -i J <fx J (03 e ^ (9 " P) j|-^(^-P) 2 ) (q,S\T°'(0)\p,S)] 

^{q,S\T^(0)\p,S) + <j>-^-{q,S\T^(0)\p,S) 



-i I d 3 q5 3 (q-p 
d(j)[(q-p) 2 



—i 



% 



Q=P 



{p,S\T°i(0)\p,S) +0(0) ^-(q,S\T°i(0)\p,S) 



Q=P 

(4.28) 



all antisymmetrized under % <-> j. 

Now 



d(f)[(q - p) 



so in the Shore- White (SW) approach 



M^{ps) = -im 



dqi 



d 



q=p 



— (q,S\T°i(0)\p,S) 



(i^J) 



q=p 



(4.29) 



(4.30) 
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which is identical to Atj^(p, 0,p',S) in Eq. (4.10), aside from normalization, which is 
irrelevant, since, at the end, it cancels out when computing actual expectation values. 

Comment 4 It is actually not possible in Eq. (4.30) to take the same S in both initial and 
final states, as in eq. (4.1). The reason is the following. A wave packet is, by definition, a 
superposition of physical states. But for a spin-1/2 particle in a physical state 

q 2 = M 2 and q ■ S = 0. (4.31) 

and we cannot integrate independently over each component of q. That is one reason why 
we have chosen to define the states using the rest-frame spin vector s instead. The correct 
procedure is then to take the same rest frame s for the initial and final states, not the same 
S. Thus the wave packet Eq. (4.25) should be modified to 

Wp,a)>= f ^-Mq ~ P) 2 U s ) (4.32) 



, (2vr)3 
and all three components of q can then be integrated over independently. 

The use of the Lorentz tensor form, as in Eq. (4.22), is central to some of the most 
interesting arguments given by Shore and White — for example their conclusion that the 
axial charge does not contribute to the angular momentum sum rules. But, as pointed out 
in Comment 3, this form is based on the incorrect expansion Eq. (4.11), and so one must 
be skeptical about their results. We shall see, however, in Section 5.2, that the conclusion 
regarding the role of the axial charge in the sum rules is correct in spite of the erroneous 
argument. 

To summarize: the Shore- White wave-packet analysis provides a justification for the 
manipulations in the Jaffe-Manohar treatment for the spatial components of the angular 
momentum tensor. But the analysis is not correct in general because it utilizes the imper- 
missible simplification S = S' and also makes an incorrect use of Lorentz covariance. 

Given these criticisms, it is interesting to consider one important practical case of a sum 
rule that can be derived from Eq. (4.22). Namely, for a longitudinally polarized proton 
moving along OZ i.e. in a state of definite helicity \ with p along OZ one has the sum rule 

\ = ±A£ + AG+ (L q ) + (L G ) (4.33) 

where AE, AG are the first moments of the polarized quark flavour-singlet and the polarized 
gluon densities, respectively, and (L q ) , (L G ) are contributions from the quark and gluon 
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internal orbital momenta. This sum rule is so manifestly intuitively right that it is almost 
impossible to contemplate it being incorrect [12]. And indeed, it is correct. This is the one 
unique case where the flaws in the deduction are irrelevant, as will be seen later. 

In Section 5 we will present what we believe to be a correct evaluation of the wave-packet 
definition of the recalcitrant forward matrix element. 

5. A DETAILED WAVE-PACKET TREATMENT 

Given our critical stance vis-a-vis the treatments of Refs. [3] and [4] it is incumbent upon 
us to proceed with caution. We shall therefore present a detailed wave-packet treatment 
for two separate situations, for a relativistic quantum-mechanical Dirac particle and for the 
field-theoretical case. It will emerge that there is complete agreement between the results 
for the two cases and that they differ, in general, from Eq. (4.23). 

We shall also address the question as to the validity of Eq. (2.12), which was assumed in 
the derivations in Refs. [3] and [4]. In order to do this we shall need to distinguish between 
-^orb ( x ) constructed out of the canonical energy- momentum density T£ (x) as in Eq. (2.4) 
and M fluX (x) used in Section 4.1 and Section 4.2 based upon the symmetrized T^ lX {x) ) as in 
Eq. (2.5). 

5-1. Relativistic quantum-mechanical Dirac particle 

We construct the wave-function corresponding to a superposition of physical states cen- 
tered around momentum p, all of which have rest-frame spin vector s: 

lM*,«) = (^3 j d ^e-^-^e^-^u( q ,s), (5.1) 



where q° = a/q 2 + M 2 , and 



"(«'*) = V -UT\ JLJL I* (S) (5 ' 2) 



9 /If I 

q° + M 

with 

X ] {s)a x(s) = s. (5.3) 
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The constant N, whose value is irrelevant for the moment, is adjusted so that ip PtS is 
normalized to 1, i.e. 

/ £x^l yS (x)^ s (x) = 1. (5.4) 

Since our aim is to provide a sensible prescription for a plane-wave state of definite mo- 
mentum p, we shall, at the end, take the limit A — ► oo. In this limit the values of q that 
contribute in the integral Eq. (5.1) are forced towards p, so we make a Taylor expansion of 
(l/g°) u(q, s) about the point p and keep only those terms that survive ultimately in the 
limit A — ► oo. 

Note that we are unable to carry out the analysis if t can be arbitrarily large. The point of 
using the wave-packet is to produce a cut-off in the divergent spatial integrals in Eq. (2.14), 
but this does not produce a cutoff in t . Thus the wave-packet approach can not be used 
if the operator densities are integrated over four-dimensional space-time. Since we are only 
interested in conserved densities, we will take advantage of the time-independence to choose 
t = in Eq. (5.1). 

It turns out to be sufficient to keep just the first two terms of the Taylor expansion, which 
yield: 



— u(q, s) = — <j u(p, s) 



P° + 2M 

2p02 u(p,s)p+ 

(q-p)} (5.5) 




p° + M ° 

\ p (p + M) 

with p° = a/p 2 + M 2 . 

The term (q — p) is removed from under the integral in Eq. (5.1) via 

(q-p)e^-^ = ^e-^-rf. (5.6) 

The remaining integral is just the Fourier transform of a Gaussian, yielding a factor 
exp(— a; 2 /(4A 2 )) exp(ip ■ x), and the V p on the RHS of Eq. (5.6) then produces a term ix. 
The result is 

i/> PtB (x, 0) = Ne- x2 / 4x2 e^ x u PjS (x) (5.7) 

where 

9 M ( 1 \ 3 , x 

N 2 = — -^=- 5.8 

p° Vv^a; v ; 
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and 





\ ix-p p° + 2M 


P,3\^J 


2A 2 2p 02 {p° + M)_ 




i lp° + M 




l 2A 2 (p° + M)V 2M 



u{p, s) + 





a ' x l^THx 



xpap I X(s)- (5.9) 

p [ \p [> + M) 



The structure of Eq. (5.9) is extremely instructive in understanding both the difference 
between matrix elements of local operators like T^ x (x) and compound ones like x p T flX (x), 
and the question as to how many terms of the Taylor expansion are necessary. 

On the one hand the Gaussian implies that the effective values of x satisfy |cc| < 2A. On 
the other, the n th term in the Taylor expansion provides a term of order (|a;|/A 2 ) n . For a 
local operator there are no other factors of x present, so even the first-order term of a Taylor 
expansion can be ignored in the limit A — ► oo. For the angular momentum, on the contrary, 
there is one explicit factor of x. The first-order term in the Taylor expansion is then essential, 
but higher order terms can be disregarded as A — ► oo. Of course, this is totally analogous to 
what happened in Section 4.1 , where the results involved a first derivative. The difference is 
that here the derivative has been calculated accurately and without recourse to an incorrect 
assumption of Lorentz covariance. 

Now that a sufficiently accurate wave function has been obtained, we turn to the cal- 
culation of the expectation values of the operators. In Dirac theory the canonical energy- 
momentum density is [13] 

Tg x (x) = -i){x)^d x i){x) + h.c. (5.10) 

and the orbital angular momentum density is 

M& x (x) = x»T£\x) - x*T&(x) (5.11) 

and the spin density is 

KtW = ^(aOyvVto- (5-12) 

The orbital angular momentum operators M^ Th and the spin angular momentum M s J p in are 
the space integrals of M or ^ and M s ^ m respectively, calculated, in our case, at t — and are 
not separately time independent. 

The calculation of the spatial derivatives of the wave function (5.7) needed in Eq. (5.10) 
and the subsequent Dirac algebra is straightforward, but laborious and will not be spelled 
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out here. Helpful is the fact that odd functions of x will vanish under integration J d 3 x since 
the non-symmetric term exp(ip • x ) cancels out in constructing Eq. (5.10). The terms in 
the spatial derivative that survive give 



&ij)j,,s{x, 0) = iNe 



- x2/4x2 e ip - x ui a (x) 



(5.13) 



where 



<«(») 



pj 



1- 



^ 



i(p° + 2M) p-x 
2\ 2 {p° + M) 2p 02 



ix- 



2X 2 < M(p ' S) 



2A 2 (p° + M) V 2M 



p° + M 





p ■ xer ■ p 



X{s)- 



(5.14) 



Then, keeping only the terms that survive in the limit A — ► oo, we find ultimately: 

1 



<b 



[p J (p x s) 1 — p\p x s) J ] 



(5.15) 



2p°(p° + M) 

or via the anologue of Eq. (2.9), we express the orbital angular momentum vector L in terms 
of the independent vectors p and s: 

1 



2p°(p° + M) 
1 



[p x (p x s) 



2 p o(pO + M ) 
For the spin angular momentum we find 



[p s-(p- s)p] 



(5.16) 



ijk 

M 13 - = — 

spin 2p0 

ijk 



Ms k + {P - 8)P 



p° + M 



= —S k M 
2p° 

where we have used Eq. (4.4). 

For the spin vector S of the system Eq. (5.17) yields 

S= V 



(5.17) 



Ms + t^P 



p^ + M^\ {5Al 

Adding Eqs (5.16) and (5.18) we have the remarkable result that the term proportional to 
p cancels out, and 



J=\s. 



(5.19) 



For later use we write this as 



J' 3 =M^ h + Mjin 



1 ijk k 
2 t a . 



(5.20) 
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In the above we have used the canonical form of the total angular momentum, where it 
is split naturally into an orbital and spin part. In Section 4, however, following Refs [3], 
[4] we utilized the angular momentum tensor built from the symmetrized Belinfante energy- 
momentum density. According to Eq. (2.12) this should yield the same forward matrix 
element as the total angular momentum in the canonical approach. For our quantum- 
mechanical example we can check whether indeed M JJI = J u . 

For the symmetrical energy- momentum density we have T^ A = \{Tq +T c m ) so that from 
Eqs. (2.5) and (5.10) 



M^\ 



x 



\M^{x) + \\ l - [x u ^{x)-f x d^{x) - x x ^j{x)Yd^{x)] + h.c. }> . (5.21) 



For M Q%: >(x) we now require the time-derivative dijj PtS (x, t)/dt\ t= o of the wave function (5.1). 
The time derivative gives a factor — iq° which is expanded as 



-m 



Po + 



p{q-p) 

Po 



(5.22) 



and via the steps explained after Eq. (5.5) this becomes a multiplicative factor: 



—i) P:S {x,t) 



t=0 



Po 



ip ■ X 

2X^0 



ip p>s (x,0). 



(5.23) 



After further laborious Dirac algebra we find 



M 13 



1 ijk k 



(5.24) 



in complete agreement with Eq. (5.20). 

Thus despite the fact that M 0l ^{x) and J°^(x) differ by the divergence of a compound 
operator, and despite the fact that the definition of the forward matrix element involves 
non- forward ones, it seems that Eq. (2.12) is indeed valid, in a wave-packet approach, as 
was claimed in Ref. [4]. 

In the next section we shall see that we obtain the same, to us surprising, result Eq.(5.20) 
in field theory, and we shall then compare it to the results in Sections 4.1 and 4.2. 

5-2. Field theoretic treatment 

Analogously to the quantum mechanical case we construct a wave-packet state \^ PtS ) as 
a linear superposition of physical (canonical spin) plane wave states of momentum q all of 
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which have the same rest-frame spin vector s: 

IT \ N 



VW) 



£ qe -^-p) 2 



|9> s ) 



(5.25) 



Note that we label the momentum eigenstates by s, not by the covariant spin vector S. As 
before we shall consider the limit A — > oo and we normalize the state to 1. 



\ ^p,s| ^ p,s/ x 



It follows that 



j\T 2 



2A 



2tt/ 2p 



(5.26) 



(5.27) 



with p = a/p 2 + M 2 . 

Let us first consider the general structure of the matrix elements of some local operator 

density 0(x) between the states \^ p , s )- We have 

N 2 r 
(^V, s |O(x,0)|* p , s > = — — Z J d 3 qd 3 q'e- x2 ^-^e- X2 ^-^ 2 e iX ^-^(q',s\O(0)\q,s) (5.28) 

where we have temporarily kept p' distinct from p as an aid to calculation. 

The matrix element on the RHS is expanded in a Taylor series in q and q' about the 
points p and p' respectively and to the required order is of the form 



i,s\O(0)\q,s) = f 1 (p,s) + (q-p) • / 2 (p, a) + (q' - p') ■ f 8 (j/,a) 



(5.29) 



Factors like q — p are transformed into -^ etc. under the integral in Eq. (5.28), so that 



(V p ,, s \O(x,0)\%, s ) 



h + h ■ ^ + f. 



V' 



2A 2 J 2A 2 



N 2 ( d 3 qd 3 q'e- x2 ^-rt 2 e- x2 «- p '^e lx <i-^ 

(5.30) 
Putting r = q — p, r' = q' — p' the integral in Eq. (5.30) becomes 



e i X .(p- P ') / d 3 re -\ 2 r* e ir. X / d 3 r , e -\ 2 r>* e -i r >. X ^g-Q 

so that V p can be replaced by ix and V p ' by — ix in Eq. (5.30). 

Carrying out the integrals in Eq. (5.30) and putting p' = p once again, we end up with 



(il> p , s \O(x,0)\il> PtS ) 
where, via Eq. (5.27) 



ix 



A(P. s ) + 777? ' (^(p, s ) ~ /a(p, *)) 



C 



2A 2 

AT 2 



c 



-x 2 /2\ 2 



1 1 1 



A 3 (x/2A) 6 (2tt)3/ 2 2p A 3 



(5.32) 



(5.33) 
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SO t licit O ~77\ — vT79" o • 

We now apply this to the case where 0(x) is the canonical energy- momentum density 
Tq V (x). The most general structure of the matrix elements of the conserved operator Tq V {jS) 

is 

(q',s\T^(0)\q,s) = u(q', «){G(gV + ?V) 

+ M(q p q" / + <A' M ) + MS[(q + q'YY + (q + q') V] 

+ (g-g'-M 2 )(G-H)(?^ 

+MA[(q + g')V - (q + q') u Y}}u(q, s) (5.34) 

where the u(q,s),u(q',s) are the usual canonical Dirac spinors normalized to mm = 1 and 
G, H, S and A are Lorentz scalars. Note that all terms, except the A-term, are symmetric in 
« <-> v. In order to identify the function f\ and the vector function f^ — fs in Eq. (5.29) 
we have to expand Eq. (5.34) about q = p and q' = p . The Gordon decomposition 

u(q')ru(q) = (l ±f^u(q')<q) + %^ u{q')a^u{q) (5.35) 

is helpful, because to the accuracy required we can replace q and q' by p in the factors 
multiplying (q' — q) p and then use 

u(p, s)a^u(p, s) = -L^ afS S aPf3 (5.36) 

where S is given by Eq. (4.4) and the convention is e 0123 = +1. Then 

(q f , s \T^(0)\q, s > = [B(q»q» + gV) + (q ■ q' - M 2 )(B - C)<T 

+C(q p q" / + q u q' p )}u(q', s )u(q, S ) 
+ ^ q '~ q)p {S[p^ puafi + pV 40 * 3 ] 
+k[p p e puaf5 - p v e p » af5 ])S a pp (5.37) 

where B = G + S and C = H + S. 
Putting 

u(q, s ) = u(p, s) + (q- p) k u k (p, s ) 
u(q', s ) = u(p, s) + (q' - p) k u k (p, s ) (5.38) 
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where u k (p, s ) = d/dp k u(p, s ) and similarly for u k , we find, after much algebra, that for 
Tq U (0) the function j\ and the k th component of fi — fs, needed in Eq. (5.32) are: 



/i = 2DpV 



(5.39) 



where 



and 



D = B + C 



(5.40) 



U'2-hf = 2Dp^[u(p,s)u k (p,s) -u k (p,s)u(p,s)} 



-2i 



(S + A) 
M 



p'' 



M F 



P Ova/3 _ kua/3 

Po 

■ fa 

P ana/3 _ kfia/3 

Po 



S a P(3 
SaP/3 



(5.41) 



so that from Eq. (5.32) 



C 



<Vw^(z,0)|Vv, s > = ^e- x2 / 2X2 {2D P y 



+i 



x 



2A 2 
2i 



2p°p i (u(p, s)u (p, s) — u (p, s)u(p, s)) 

k 

M r >' 



(S_^) /( ^ e o^_^ )5 ^ + 2i (S-A).., 



M 



pPe k0 ^S a p p 



(5.42) 



Consider now the matrix element t 



d 3 xx l {^ s \T^(x,0)\^ s ) - (i^j) (5.43) 

We see that integrating over space kills the first term in Eq. (5.42), and the second only 
contributes if k — %. The spatial integral is then simply 



|/^.)v— = £ 



(5.44) 



t Note that here, in the general field theoretic case, M or b is actually the sum of the quark orbital angular 
momentum plus the full angular momentum of the gluons 



26 



where we have used Eq. (5.33) for C. The A 2 in Eq. (5.44) cancels the remaining 1/A 2 in 
Eq. (3.42), so that we can take the limit A — *■ oo and have 

-u l (p, s)u(p, s)} - ^ [p l e ^ - p e^ + p j e 0i ^]S a p f3 

j IK 

Part of the term multiplying S is symmetric under (i *-* j), so cancels out. The other terms 
multiplying S and A are antisymmetric so just get doubled under (i <-> j). Also 

u(p, s)u l (p, s) — u l (p, s)u(p, s) 

= u(p, s)u l (p, s) - [u(p, s)u l (p, s)]* 

= 2ilm[u(p, s)u l (p,s)] (5.46) 

Thus eq. (5.45) becomes 

{i> P ,s\ M oJi>P,s) = -D[p>7m[w(p, s)u l (p, s)} - (i <-► j)} 

—J^SaPp + ■^-(P i z 0iafi - P 1 ^ - Po^)S a p? (5.47) 
Finally we expand the spinors as in eq. (5.38) and find that 

Im[u(p, s)u\p, s)] = -—^——^ lm p lSm (5.48) 

and comparing Eq. (4.13) with A o (0) = 2 [(see Eq. (4.16)] with Eq. (5.37), we see that 
B + C = 1 and so, via Eq. (5.40), D = 1. Thus 

& P ,s\MM,s) = 2M(^TM) [P ' (PXS)J ~ pJ(PXS)i] 

-j^e^SaPp + ^(p** ** ~ V 3 ^ - Poe ija ")S a pp (5.49) 

We shall return to consider the interpretation of this result presently. First, though, we 
use it to deduce the structure of the matrix elements of J l \ which, by Eq. (2.12), are the 
same as those of M l K The latter are built from the symmetric energy momentum density 
T^ v = \{Tq / + Tq 1 ). Since in Eq. (5.34) all terms are symmetric in //, v except the A-term 
which is antisymmetric, we obtain the matrix elements of T^ v by simply putting A = in 
Eq. (5.49). Thus 

<^, s |M^ p , s > = 2M{p l +M) [P l (P x sy - pP{p x «)'] - ^e^S a pp. (5.50) 
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We obtain the value of S by choosing a wave-packet with p = (0, 0,p) and s = (0, 0, 1). 
This is then a helicity state lV'1/2) and should be an eigenstate of J z with eigenvalue 1/2. 
Hence for this state 

ftM-Wi/2> = (Vi/ 2 |M 12 |^ 1/2 ) = 1/2 (5.51) 

which implies 

\ = "-|e 1203 (<SoP3 - <%>o) = j^(SoP 3 - S 3 Po ) (5.52) 

Now from Eq. (4.4), for the present case, 



S = p/M 
S 3 = 1 + 



P 



Hence Eq. (5.52) becomes, 



M(p + M) 



\ = w^ p2 - pl) 



Po/M 



Putting this into Eq. (5.50) gives 



1 



2M {po + M 

+6 iJa ?S aPp 



\p\px 



\3 _ n) 



p J (p X s 



(5.53) 



(5.54) 



(5.55) 



Let us first compare this result to what we obtained in Eq. (5.20), Eq. (5.24) for the 
relativistic quantum-mechanical Dirac particle. The last term in eq. (5.55) can be written 



e ijap S a 



P0 



c im 



4jk 



c ijk 



c ijk 



[Sopfj - Sppo] 

}" sk ~ P ^ h 



Pqs 



M(p + M) 



P' s k 

■P 



P 



P ' s k 

- V 



Po + M- 



Putting this into Eq. (5.55) we have, finally. 



{^ S \M^\^ S ) 



2M [po + M 



p X (p X s) 



} k +Pos k 



P * s k 

— V 

Po + M F 



e ijk 

2M 

e ijk 



Po 



p + M 



(5.56) 



(5.57) 
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This is exactly the result found in Eq. (5.20), Eq. (5.24). Taken in conjunction with 
Eqs. (2.9), (2.10) and (2.12), the result Eq. (5.57) corresponds to the first term in Eq. (1.1). 
The derivative term in the latter equation is missing here because of our use of a wave 
packet. We shall compare Eq. (5.57) to the results of [3] and [4] in Section 6. 

Let us return now to the canonical form of the angular momentum. It is shown in [4] for 
QCD that aside from the terms which vanish by virtue of the equations of motion or which 
give no contribution to forward matrix elements 

M^\x) = M^\x) + M%*(x) (5.58) 

where in the convention e 0123 = +1, the axial density is 

M^(x) = \e^ x ^{x) lbla ^x) (5.59) 

and a sum over the color labels of the quark fields is implied. 

Note that Eq. (5.59) is not the same as the spin density in Eq. (5.12), for arbitrary values 
of /i, z/, A. But Eq. (5.12) and Eq. (5.59) are equal if \x ^ v, A. Thus 

M^Jx) = M^ix) (5.60) 

where "spin /" means the fermionic spin terms in the angular momentum operator. This is 
the basis for the result in Eq. (2.12). The axial density is a local operator, so that 

(p',s\Jd*xM^(x)\p,s) = (2nf5"(p'-p)(p',s\M^(0)\p,s) (5.61) 

The only possible structure for the matrix element on the RHS is 

(p,s|M^(0)b, S > = -aoe^ (5.62) 



where ao is known as the axial charge. Thus for the expectation value 

(p',s\j<PxM^ l {x ) Q)\p ) s) = _ao_ eflvXps 
(p', s\p, s) 2p p 



(5.63) 



From Eq. (5.60) we thus have for the spin density 

n "-e ijk S k (5.64) 



(p',s\fd 3 xM°;{ nf (x,0)\p,s) __ „ 



(p 1 , s\p, s) 2p 
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After a little algebra, using Eq. (4.4), we write the fermionic part of the spin vector S^ of 
the nucleon in terms of the independent vectors and p and s:' 



of _. °o 
2p 



Msk + ^-J7P k 



(5.65) 



Similarly, for the vector representing the fermionic orbital angular momentum vector plus 
the full gluon angular momentum, Eq. (5.49) yields 



1 , ,c (1 >M\_ , p-s 



L ' + *=I2-VJ S '- A ^TM) P ' (5 ' 66) 



Adding Eq. (5.65) and eq. (5.66) gives 

3 
Then Eq. (5.57) implies that 
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M , 

-- A- 
Vo 


a 
~~2~ _ 


s - (A- 


a p-s 
2 J p (po + M) 


lat 




* = ? 





p (5.67) 



(5.66 



and the contribution from the axial charge cancels against the antisymmetric term in the 
fermionic orbital plus full gluonic angular momentum. 

This supports the claim in [4] that the axial contribution totally cancels out when taking 
forward matrix elements of Eq. (5.58), at least for the Oij elements of M^ vX . So although the 
proof given in [4] is invalid because it suffers from the problems mentioned in Comments 1 
and 3, it is correct in the contention that the axial anomaly cancels out of the angular 
momentum sum rules. ' As the authors of [4] realize (cf. their Eq. (21)), this must oc- 
cur: because of the non-renormalization of the energy-momentum tensor, the total angular 
momentum must not be anomalous. See also [14]. 

6. AN INDEPENDENT APPROACH 

We have argued that the forms for the angular momentum tensor given in [3] and [4] are 
incorrect, in spite of their appearance of explicit Lorentz invariance. On the other hand our 



t Note that we are deliberately not calling S-* the quark spin contribution to the nucleon spin. In an 
interacting theory with an anomaly the connection between S^ and the spin carried by the quarks is 
subtle and will be discussed in Section 7. 

' This does not mean that the spin carried by the quarks does not contribute to the spin of the nucleon. 
As explained in [8] ao contains also an anomalous gluon contribution. 
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procedure has led to the surprising result Eq. (5.57) that for a spin | particle described by 
a canonical spin state 

(p', s\M ij \p, s) = 2p (e ljk s k /2 + ip % & - ijPd 1 ) (2tt) 3 5 (3 V - p), (6.1) 

though the derivatives of the delta-functions shown here disappeared in our wave-function 
treatment, via partial integration, for the kind of wave packet chosen in Section 5. The 
key part of Eq. (6.1), namely the first term on the RHS, was derived in several different 
ways, where always careful attention was paid to the definition of matrix elements required 
by the compound nature of the angular momentum density operator. However, it does not 
share the explicit Lorentz invariance of the form Eq. (4.23) and so one may question its 
correctness. 

In this Section we will confront this issue by deriving Eq. (6.1) by a totally different 
method, valid for arbitrary spin, which is based on the rotational properties of states, and 
which circumvents completely the use of the energy-momentum tensor, and demonstrate 
that Eq. (6.1) is exactly what is expected on very general grounds. We shall then show that 
our result does have the correct Lorentz transformation properties. Finally we shall derive 
the analogue of Eq. (6.1) for particles described by Jacob- Wick helicity states. 

6-1. Canonical spin state matrix elements 

In order to utilize the rotational properties of the canonical or boost spin states we need 
to display explicitly the Wigner boost operators used in defining the states of a moving 
particle in terms of the rest frame spin states quantized in the z-direction |0, m). Then the 
transformation properties of the states becomes explicit. 

We will use the definitions of the Lorentz group generators given in Weinberg [15]. In 
this section we will write the 3-vector operators in terms of the integrated tensor M l i 

J t = y ijk M jk (6.2) 

Kt = M i0 . (6.3) 

where on the left hand side % — x,y, z. For a particle of mass M in motion the boost (or 
canonical) state is defined by 

\p,m) = B(v)\0,m) = exp(i(p ■ K)\0,m), (6.4) 
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where v = p/po, cosh£ = po/M, and p is the unit vector along p. 

Now consider a rotation about axis-i through an angle f3. The unitary operator which 
effects this is given in terms of the angular momentum operator Jf 

Ri(P) = exp(-i/3J i ) (6.5) 

and for a particle of arbitrary spin-s 

Ri(P))\P,rn) = \RmP,n)V: m (R w (p,[3)). (6.6) 

where R\y(p, P) is the Wigner rotation. In a slight abuse of notation, we will use the symbol 
R to denote both the unitary operator for rotations in Hilbert space and the corresponding 
rotation matrix in Minkowski space. The same will be done for boosts, B. For a pure 
rotation the Wigner rotation R\y is very simple 

R w (p,f3)=R i (P); 

independent of p ; cf. for example [8]. Therefore, 

(p',m'\R t ([3)\P,m) = (p / ,m / | J R i (/3)p,n)P n s m (i?,(/3)) (6.7) 

= 2p (27r) 3 5 (3 \p' - RitfMV^RiiP)), 

using the conventional normalization 

(p',m'\p,m) = 2 Po (27i) 3 5^\p' - p)5 m , m . (6.8) 

Thus 

d 
(p',m'\Ji\p,m) = i—(p',m'\(Ri(l3)\p,m)\p =0 (6.9) 



= 2po(2tt) 
Now [9] 



d „ .d 



KijkPjQ^Sm'm + i-QpVm'miWm^o 



5^(p' -pX-6-10) 



d_ 



i— V^ m (Ri(P))\ = (Si) mln (6.11) 



where the three (2s + 1) dimensional matrices S$ are the spin matrices for spin-s which 
satisfy 

[Sj,S fc ] = ie jk iSi (6.12) 
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Thus, our final result for the matrix elements of the angular momentum operators for arbi- 
trary spin, from Eq. (6.10), becomes 

d 



(p',m'\J i \p,m) = 2p (2vr) 



&i + 1-tijkPj ' 



5 {3 \p'-p). (6.13) 



dp k _ 

For spin |, of course, the Sj are just \ times the Pauli matrices <Tj. For arbitrary spin they 
are still very simple: 

1 
(pxjm'm — n |p ( s ; m ) "m',m+l + ^ ( s i ~~ m ) "m',m-l\ 

— 1 
(Sy)m>m = "Z" [ C (s, m) 5 m ', m+ l - C(s, -m) 5 m ',m-l] (6.14) 

where 



C(s,m) = -\/(s — m)(s + m + 1) (6.15) 

For the case of spin | Eq. (6.13) is exactly equivalent to the result quoted in Eq. (1.4) and 
which we obtained in Section 5 after much labor using the wave packet approach. It is 
completely general. The second term will vanish if integrated over symmetric wave packets. 
However it must be kept for analyzing the transformation properties, as we will see, and 
must, as usual, always be interpreted in the sense of partial integration. It is very easy to 
verify that the form Eq.(6.13) satisfies the usual commutation relation relations and so is 
consistent with rotational invariance. 

Combining the result Eq.(6.13) for the case of spin \ with Eq. (1.3) leads directly to the 
result quoted in Eq. (1.1). 

6-2. Lorentz invariance 

We shall now demonstrate that, despite appearances to the contrary, Eq (6.13) is consis- 
tent with Lorentz invariance. Because of the complicated algebra involved, we shall present 
the proof just for spin | of mass M. Under Lorentz transformations, M t0 = —M° l = Ki are 
brought in so we need the matrix elements of K^. These can be obtained just as those for 
the angular momentum. For a boost of magnitude and direction oj 

S(w)=exp(iwlf). (6.16) 
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Here we are using a natural shorthand for a boost Eq. (6.4) with velocity v = uj tanh(cj) rs uj 
for small |w|. Proceeding as before we have 

(p',m'\B(uj)\p,m) = (j/,m'\B(«>)p,n)Vl£(R w (p,w)). (6.17) 

The Wigner rotation for this case is defined by 

B(B(uj)p)R w (p, uj) = B{uj)B(p) (6.18) 

and is much more complicated than for rotations. However, we will only need the Wigner 
angle for small Lorentz transformations for our discussion. By doing the matrix multiplica- 
tion explicitly for small w we find the Wigner rotation is given, in magnitude and direction, 
by 

BXW . 

a = — — . 6.19 

Po + M K J 

In this case the Wigner rotation angle depends on p and this leads to the necessarily 
more complicated commutation properties involving the boost generators. Differentiating 
Eq.(6.17) with respect to uj and then putting uj = produces 

1 (p x <x) 



{p,m'\Ki\p,m) = 2po(2tt) 



-ip di 



5 {3) (p'-p). (6.20) 



2 po + M 

We will now verify that the matrix elements of Ji and Ki transform among each other 
correctly. The Lorentz transformation of Ji under a boost uj requires that 

{p',m'\Ji\p,m) = (B^p^n'lB^JiBiujy^B^p^n) 

VlZ,(R w (p',uj))V^ m (R w (p,uj)). (6.21) 

Notice that the two Wigner rotations appearing in Eq.(6.21) are not the same as long as 
p t^ p'. It is essential to keep this is mind because of the derivatives of S- functions that 
enter. If we write this out to terms linear in uj and use 

[Ji, Kj] = ie ljk K k (6.22) 

Eq.(6.21) becomes 

(p > ,m > \J i \p,m) = (B(uj)p' ',n'\Ji + ei jk UjK k \B(uj)p,n) 

+ 1 - i^e abc a a — — 1 + i^e abc a a — — . (6.23) 

V Po + Mj m , n , V Po + Mj nm 
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We must also expand the matrix element to first order in uj. Using Eq.(6.20) and Eq. (6.13) 
with Sj replaced by cTj/2, and recalling that poS^ (p 1 — p) = p'^^^p' — p) is invariant under 
Lorentz transformations we find that 



(B(u>)p',n'\Ji + e ijk UjK k \B{uj)p,n) 



ic 



U 



k Pjd k + \oi + \ 



1 (w x (p x a)\ 



Po + M 



2 Po {2tt) 3 5^{p' 



Ph 



(6.24) 



If one inserts Eq.(6.24) into Eq.(6.23) and evaluates the terms proportional to uo one easily 
finds three components: one coming from the second term in Eq.(6.24) combined with the 
rotation matrices in Eq.(6.23), one coming from the third term in Eq.(6.24) , and one from 
the first term in Eq.(6.24) acting on the rotation matrices in Eq.(6.23). It is essential to 
carry out the derivatives in this last component before setting p' = p: 

(wx(pX(7))j .[«r-(wxp),(7i] 



2(p + M) 



+ie ijk p j d k 



4(p + M) 
(uj x p' 



1 



(T ■ (<jJ X p) 

2{po + M) 



(6.25) 



2(p + M) 

Use of the commutation relations of the Pauli matrices and the vector double cross product 
identity shows that the sum of these three pieces vanishes leaving just the matrix element 
of Ji as was to be shown in order to satisfy Eq.(6.21) or (6.23). 

Using the same techniques, one can show that the matrix elements of Ki as given in 
Eq.(6.20) also transform correctly under boosts and so Eq.(6.13) and Eq.(6.20) form a rep- 
resentation of the Lorentz group. 



6-3. Helicity state matrix elements 



We now turn to the case of helicity states which have some rather surprising properties. 
One can proceed just as here; the main difference is that the Wigner rotation becomes a 
Wick helicity rotation, always about the z-axis. This simplifies things somewhat; all the 
complication is in calculating the angle that results, the analog of Eq.(6.19). The result 
is also convention dependent, depending on whether one uses the original Jacob and Wick 
definition [7] or the later one due to Wick [10] [see Eqs. (3.2) and (3.3)] . We give here 
the result for the first case. The result of this messy calculation is that, for p = (p, 8, <f) 



(p', A'| J t \p, \) JW = (2n) 3 2p [\ Vt + i{p x V),] 5 (3 \p' - p)5 y> 



(6.26) 
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where 

rj x = cos(0)tan(6»/2), r) y = sin(0) tan(0/2), 77* = 1. (6.27) 

Although these components look a little odd — the singularity at 9 = i\ results from the 
ambiguity of Jacob and Wick helicity states at that point — it is easy to verify some important 
properties : they are manifestly diagonal in A, which is required since rotations preserve the 
helicity, and they satisfy 

(p', X'\p ■ J\p, \} JW = X 2p (2tt) 3 5 (3 V - p) 5 yx (6.28) 

and no orbital angular momentum piece survives as expected. 

It is enlightening to consider these amplitudes from a different direction: comparing the 
definitions of canonical (boost) states to helicity states we have for the case of spin | 

\p,s) = \p,m)V l J* /2 {R{s)) 

V2/ D -i/ \\<r>i/2 



\pA)jwV^(R-\ P ))V^ /2 (R(s)) 

,1/2/ 

Al/2^ 



\pA)jwVl[ 2 JR-\p)R{s)). (6.29) 



This has the appearence of an ordinary unitary change of basis, but because of the compound 
nature of J, when we apply this to the canonical form, using the spin | version of Eq. (6.13), 
we get 



1/2 / D / vv^l/2 



(p',\'\J t \p,\) JW = (2nY2p V^(R(p))V^ x ,(R(p')) 



ieijkPjdk + - Oi 



8^{p'-p). 

m 

(6.30) 

We cannot use the unitarity of the X>'s because p ^ p', and we must first pass the first 
T> l l 2 (R(p)) through the derivative before setting them equal. This produces an extra term 

-(2n) 3 2p te ljk V^l f (R(p')r Pj d k V^(R(p))5^(p'-p) (6.31) 

which is tedious to evaluate in the general case. The result of this labor is identical to 
Eqs.(6.26, 6.27). 

7. APPLICATIONS : SUM RULES 

Equipped with the expressions for the matrix elements of Ji derived in Sections 5 and 
6 and summarised explicitly in Eqs. (1.4), (1.5) and (1.6), we will derive the general form 
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for angular momentum sum rules for the nucleon and, in particular, will derive a new sum 
rule for transverse polarization. This differs from the sum rule that would be derived from 
Eq. (4.23). 

7-1. The matrix elements of J: the J z sum rule 

The first term in our result Eq. (1.1) differs from the results of Jaffe and Manohar [3]. If 
we rewrite their expression Eq. (4.23) in terms of the independent vectors p and s, we find, 
for the expectation value 



«>- = iw i (3p ° - M)s ' - 7^w {p ' s)p 7 (7 ' 1) 

to be compared to 

(Ji) = fa (7.2) 

arising from the first term in Eq. (1.1). In general these are different. However, one may 
easily check that if s = p the Jaffe-Manohar value agrees with Eq.(7.2), while if s _L p they 
are not the same. 

The agreement for s = p is consistent with the much used and intuitive sum rule 

\ = IAS + AG + (L q ) + (L G ) (7.3) 

based on the matrix elements of J z , for a proton moving along the z- axis with helicity 
A = 1 (this, as explained in Section 3, coincides with a canonical spin state), relating the 
component along p of the spin and orbital angular momentum of the quarks and gluons to 
the helicity of the proton. 

7-2. General structure of sum rules 

Consider a nucleon with momentum along OZ, p = (0,0, p), in a canonical spin state 
with rest-frame spin eigenvector along s, where s could be longitudinal sl or transverse St- 
Sum rules can be constructed by equating the expression Eq. (1.4) for the nucleon matrix 
elements (p',m'\Ji\p,m) with the expression obtained when the nucleon state is expressed 
in terms of the wave functions of its constituents (partons; quarks and gluons). 
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There is great interest in sum rules in which the partonic quantities can be related to 
other physically measurable quantities. The classic example of this is Eq. (7.3). We will 
now investigate other similar possibilities, using Eq. (1.4) as the relevant starting point. 

We have stressed the importance of a wave-packet approach in order to deal with the 
derivative of the delta-function in the equations above . As it happens, however, when 
constructing sum rules, the expression in terms of constituents automatically produces a 
term which cancels the delta-function, irrespective of the actual model wave-functions used. 

The nucleon state is expanded as a superposition of n-parton Fock states, wherein, for 
the purpose of showing the structure of the sum rules, we will not display flavour and colour 
labels. We use the "Instant" form rather than the commonly used "Light-Cone" form [16] 
since it is more suitable for discussing rotational properties. We will use the original p — ► oo 
limit in order to obtain the parton model sum rules. 



\p,m) = [(27r) 3 2p ] 1/2 ^^| 



d 3 k x d 3 k n 



xi/j piTn (k 1 ,a 1 ,...k n ,cr n )8 ( - 3 ' ) (p-k 1 -... -fc n )|fc 1 ,o-i,...fc n ,o- n ). (7.4) 

where <7j denotes either the spin projection on the z-axis or the helicity, as appropriate. ipp, m 
is the partonic wave function of the nucleon normalized so that 

J2 / d 3 k ± . . . d 3 fc n |VVn( fe i> ff i> - fe n, Vn)\ 2 S (3) (p -k 1 -...-k n ) = V n . (7.5) 

with V n denoting the probability of the n-parton state. 
The n-parton contribution is then 

fr/, m'\Ji\p, m)™ ton = (27r) 3 2p ]T J[d%] . . . [d%] [d 3 h] . . . [d 3 k n ] 



^*p,m( k l, °"l> ••• fc n> °n)( fc l> °1> -, K, °n\ J i\ k U 01, -, fen, <T n ) (7.6) 

ip pm {k 1 ,a 1 ,...k n ,a n )5 {3 \p' - k[ - ... - k' n )5 {3) (p - k ± - ... - k n ) 

where we use the notation 

d 3 k 

[ J ~ v/(27r) 3 2A; ' 

We take for the Fock-state matrix elements 

(fei, a[, ..., k' n , a' n \ J t \k u a,, ..., k n , a n ) = J^iK, <\MK, <?r) H(2ir) 3 2k°5^ {k[ - h)5 a[ai . 

r Ijtr 

(7.8) 
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so 



p',m'\Ji\p,m) = {2tc) 3 2 Po J2J2J2 f \^K\ d 3 h. . .[d 3 k r ] . . .d 3 k n (k' r ,a' r \Ji\k r ,a r ) 



n,r Hi u 



rp',m'l" ; l! ^li ---"V' ®ri •••^ni (T n)' , Pp,m{f i li G\i ■■■K r: <J r , ...fc n , <J n ) 

$m( p > _ fel - fe 2 ... - fe;... - k n )5 {3 \ P - fex - fc 2 ... - fc r ... - fe n ). (7.9) 



After some manipulation this can be written as: 



p'XlJ^m) = (27r) 3 2p ^^ f d 3 kd 3 k'5 (3 \p , -p + k-k')p^^(k , ,k) a 



n a,a 



, l (k'.a'UAk.a)—, l (7.10) 

where we have introduced a density matrix for the internal motion of type 'a' partons in a 
proton of momentum p : 

P "j a m (k',kr = £EE<w<w 

/ d 3 k' r d 3 k l ...d 3 k r ...d 3 k n 5 {3 \k - k r ) 5^ 3 \k' - k' r ) 

V'p'm'l"'!' °"l' "•">' °V> ■■■k n , (7 n ) 1ppm{kl, 0"l, ...K r , 0" r , ...K„, <7 ra ) 
<J (3) (p-fel...-fe r ...-fe„) (7.11) 

Here a, which we will frequently suppress, denotes the type of parton: quark, anti-quark 
or gluon. The sum goes over all Fock states and, within these states, over the spin and 
momentum labels r corresponding to the parton type a. Eqs. (7.10) and (7.11) are the basis 
for the angular momentum sum rules. 

The two terms in Eq. (1.4) applied to the parton matrix elements in Eq. (7.10) suggest a 
spin part and an orbital part for quarks and gluons. This decomposition is a little misleading 
at this level, as we will see, but we will use it here to organize the various pieces. First 
consider the spin part of the matrix element when k is the momentum carried by a quark. 

(p f , m'\Ji\p, m ) quarkspin = (2tt) 3 2p 5 (3 \p' - p) j d 3 kd 3 k' <5 (3) (fc - fe') 

xEiW^**'.*)'. ( 7 - 12 ) 

<t,<t' 

where here ctj denotes the Pauli spin matrix of Eq. (1.4). 
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The spin part for the gluons is completely analogous, but now a and a' in Eq. (7.10) refer 
to the gluon helicity A. From Eq. (1.5), which is diagonal in helicity, we obtain 

(p', m'\Ji\p, m ) gluonspin = (2tt) 3 2 Po 5 {3) (p' - p) I d 3 kd 3 k' 5 (3 \k - k') 

THX ptx m {k\kf (7.13) 

The orbital part is somewhat different because of the derivative of the 5-function that 
enters, and we have stressed the need for a proper wave packet treatment, as carried out in 
Section 5. However in Eq. (7.9) or (7.10) the partons are not in plane wave states and the 
partonic wave function ip plays the role of a wave packet. Thus we may proceed directly by 
inserting the orbital piece of Eq. (1.4) into Eq. (7.9). 

We have 

(p',m'|J«|p,m> arbital = (2n) 3 2 Po J2J2Yl I ' ^K^h . . .d 3 k r . . .d 3 k n 

n { CT } r 
' l Pp',m'\kli °"l) •••> "V> ®ri ■■■fcn, O'n) Yp,m\kl, 01, ..., fc r , <7 r , ...fc n , <J n ) 

SWtf -k x - ... -K- ...k n )S^(p - fex - ... -K- ...k n ) 
i{k r xV kr )^ 3 \k' r -k r ) (7.14) 

Integrating over k r by parts yields 
(p',m'\J t \p,mr h ^ 1 = -(2n) 3 2p J2J2J2 [ d%d% . . .d 3 k r . . .d 3 kj ^. 

n {a} r J V K r 

rp',m'("'l> °"l5 •••; "V> ®ri •••" ; n? ®n) 

5 {3 \p' - fex... - k' r ... - k n )5 {3 \k' r - kr) (7.15) 

lyK r X V kr)i[Yp,m{kl, c r .., K r , & r ..k n , u n )o [jp K\.. K r .. K n )\ 

The derivative produces two terms. The one arising from the derivative of the delta- 
function is 




= -(27r) 3 2p ^^ / d 3 ki d 3 k n i)* pml (ki, a 1} ..., k n , a n ) ip p , m (ki, 01, ..., k n , a n ) 

n {a} J 

S(V(p> - kl ... - k n )J2^K x V fer )^(p - fei... - fen) (7.16) 

r 

Now it is easy to check that 

J^(fe r X V fcr )^ (3) (p - fei - .- - fer - ... - fen) = ~{p X V p )^ (3) (p " fei ~ .- - fen) (7.17) 
r 
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and putting this into Eq. (7.16) and using the normalization and orthogonality of the wave 
functions, this term produces 

d 
2p (27r) 3 te ljkPj —5 i3 \p' - p) b mm , (7.18) 

which just cancels the derivative of the delta-function in Eq. (1.4). 
The other term in the differentiation in Eq. (7.16) yields 

2p Q {2-K)H^{p r -p){L l ) a m , m (7.19) 

where (L i )^ n , m is the contribution from the internal angular momentum arising from partons 
of type a, given by 

( L i}m'm = Yl Yl / d3kl d3k n%,m'(kl, <?1, •••, K, CT n ) 

n {a} J 

y^{[-j(fc r X Vfe r )i]V'p,m(fel,0 : i,...fe r ,(T r ,...fe n ,Cr n )} 

r(a) 

5^(p-k 1 -...-k n ) (7.20) 

where the sum over r(a) means a sum over those r- values corresponding to partons of type a. 
Note that a can refer to both quarks and gluons; the structure of Eq. (7.20) is the same for 
both. Note also that the orbital angular momentum defined in this way is not in general the 
same as that given by the matrix element M^ b , Eqs. (5.49) and (5.66), which contain gluon 
spin-dependent parts which here are included in the spin part of the matrix elements of J,. 
This difference is important for transverse polarization but not for longitudinal polarization. 
Furthermore, it is important to realize that the orbital angular momentum defined in this way 
depends on the basis states used for the partons. In particular, because of the momentum 
dependence of the transformation from canonical to helicity basis, Eq. (6.29). 

p™' A m (fc / , k) = V%{R{k')Y p^r mr {k\ k) Vl' 2 rX {R{k)) (7.21) 

the orbital angular momentum will not be the same in the two bases; cf. Eqs. (6.30) and 
(6.31). Of course, for k z — > oo with finite k x , k y the z-component will be the same but the 
other components will not be. 

In this discussion we have used a fixed-axis quantization for the quarks since it has a 
more transparent connection to the polarization states of the proton, which in most cases 
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are more naturally described in that way. At the same time, because the gluons are massless 
it is most natural to use helicity states to describe them, so we are led to a mixed notation. 
There is no real problem with this, but there may be occasions where one wishes to treat 
the quarks in helicity states as well. 

Putting Eqs. (7.19), (7.12) and (7.13) into Eq. (7.10) , utilizing Eq. (1.4) for its LHS , 
and cancelling the factors 2po(2n) 3 5(p / — p), we end up with the general sum rule for a spin 



\ nucleon: 




2\&i)m'm = J CI K 


"SfoW ftf?(k,k)*« + \ mik) pf x m (k,kf 


'\^i)m'm ' {■'-'i/m'm 


where rji is given in Eq. (6.27). 





(7.22) 



7-3. A new sum rule 

For proton matrix elements of J z , Eq. (7.22) is non-vanishing only for w! = m. It then 
becomes the classic sum rule Eq. (7.3). There is one other independent sum rule that can 
be obtained from this general one; one way to obtain it is to consider the matrix elements 
of J x which are non-vanishing only for w! = —m. The LHS of Eq. (7.22) is then equal to \. 
The quark spin contribution to the RHS (an identical expression holds for the antiquarks) 
is 

\ J ' Sk\[p + + -_ + pt+ + pl+ + p- + t] q (7.23) 

where +/- refers to ±|. By rotating the system through n about the z-axis, it is easy to 
see that elements of p™, '™ with (— l) m ~ m -<r+<* — _1 are Q ^ under this rotation and so will 
integrate to zero when integrated over &t- This enables us to rewrite the expression (7.23), 
the quark contribution, in a way that has a nice interpretation, viz. 

\J ' d z k\[pX+_ + p + _X + p- + Z + pZ- + + pXz + pt- + + pZX + p- + t] q . (7.24) 

Consider the proton state with spin oriented along OX, perpendicular to the direction of 
motion 

\p,a x ) = A={\P,m= 1/2) + \p,m = -1/2)} (7.25) 

To understand the content of expression (7.24) write schematically 

/#'" = Yl O', X) ip m (cr, X) (7.26) 

X=all 
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Now the number density of quarks with spin along or opposite to OX, denoted by ±s x in 
a proton spinning along OX is 

q±s x /sM= E \^s x (±s x ,X)\ 2 (7.27) 

X=all 

where 

lM±*») = i[^ + (+) ± M-) + V-(+) ± V-(-)] (7-28) 

so that (suppressing the ^x=a«) 

fe/ s ,(fc)-^V^( fc ) = ^{[^ + (+)+^-(+)]*[^(-)+^ + (-)]} (7-29) 

which, via Eq. (7.26), is exactly the integrand in Eq. (7.24). Thus the expression (7.23) is 
equal to 

\\ d 3 k [ q s - x/Sx (k) - q-- Sx/Sx {k)\ = \ J dx d 2 k T [ q s ~ x/Sx (x, k T ) - qs x / Sx {x, k T ) ] (7.30) 

and there is an an analogous term for the antiquarks. 
Now it is known [19] that 

ql /Sx (x,k T )-qU x/Sx {x,k T ) = A' T q a (x,k 2 T ) 



+ cos2(j)^h^(x,k 2 T ) +sm<j)j^hi a (x,k 2 T ) 



(7-31) 

where is the azimuthal angle of k? , A' T q a (x, k^) is the same as h®(x, k 2 -) in the notation 
of Ref. [17], and 

A T q a {x) = h1(x) = I d 2 k T A' T q a (x, k 2 T ). (7.32) 

Substituting Eq. (7.31) into Eq. (7.30) and integrating over the direction of kx , we end 
up with the quark contribution to the RHS of Eq. (7.22): 

\ I dxh±{x) = \ I dxJ2 ha i(x) = \ I ' dx^A T q a {x). (7.33) 

a,a " a,a 

We turn now to the gluon contribution to the RHS of Eq. (7.22), which is 

d 3 k %(*;) (i [p+ + - ptt-i + Pii ~ Pir-i + Pii~ - P-i-i + Pii - P-t-i ]), (7-34) 



where ±1 refers to the gluon helicity. Once again we have added in terms which integrate 
to zero in order to get a nice interpretation in terms of densities. (Recall that i] x contains 
the factor cos</>, and the factor p™ '™ with (— l) m _m_A + A = +1 are even under <fi — ► it ± 0.) 
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Now consider 

AG h/Sx = G 1/Sx - G-!/ Sx 

= £{|^(i,*)l 2 -liU-i,*)l 2 } ( 7 - 35 ) 

X=all 

Carrying out the analogue of Eq. (7.28) we find that the RHS of Eq. (7.35) is exactly equal to 
the terms in parenthesis in Eq. (7.34). Thus the gluon contribution to the RHS of Eq. (7.22) 
is 

r d 3 k Vx (k)AG h/Sx (k) = fdxd 2 k T Vx(k)AG h/Sx (x, k T ) (7.36) 

It is easy to see, geometrically, that AGh/ 3x (x, kx) contains a factor k x and we make this 
explicit by writing [17] 

&G h/Sx (x, k T ) = ^9?t(x, k 2 T ). (7.37) 

Then the contribution of the gluon spin to the RHS of Eq. (7.22) i.e. to the proton whose 
spin is in the x-direction is 



AG h/Sx = dxd 2 k T r] x j^g^ T (x,k] 



vr j dx k T dk T ^ XP ±p XP g? T (^ 4))- (7-38) 



Yp, m > Yp, s x /-- 



^ p l +^„_1 



(7.39) 



where we have used Eq. (6.27). As p — ► oo this piece vanishes and so the gluon spin does 
not contribute to the transverse spin sum rule. 

Finally, the internal orbital angular momentum terms {L x )\ and (L x )^ are obtained 
from Eq. (7.20) by the replacement 

1 

2 v,-2 

Putting together the various pieces of the RHS of Eq. (7.22) we obtain a new, transverse 
spin sum rule. Since the same result holds when considering J y with the proton polarized 
along OY, we prefer to state the result in the more general form: for a proton in an eigenstate 
of transverse spin with eigenvector along St 

3 = 3 E JdxA T q a (x) + E (Ls T ) a (7.40) 

<?>5 q,q,G 

where L ST is the component of L along St- No such sum rule is possible with the Jaffe- 
Manohar formula because, as p — > oo, Eq. (7.1) for i — x,y diverges. 
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This has a very intuitive appearance, very similar to Eq. (7.3), but there are a couple of 
points to be made. First of all, this sum rule is distinct from the tensor charge 5 sum rule 
of Jaffe and Ji [18] because here the quarks and anti-quarks add whereas in S the difference 
enters. It is possible that this difference may be of some use in disentangling the quark and 
anti-quark transverse spin structure functions. Second, the integral does not correspond to 
the matrix element of the axial vector current, which via Eq. (5.65), vanishes for p — ► oo 
when p ■ s — 0. 

The structure functions A^g a (x) = h\(x) are most directly measured in doubly polarized 
Drell-Yan reactions where the asymmetry is proportional to 

J2el[A T q a (x l )A T q a (x 2 ) + (1 <-> 2)]. (7.41) 

a 

For a detailed discussion see Ref. [19]. 

They can also be determined from the asymmetry in semi-inclusive hadronic interactions 
like 

p + p(s T ) -> H + X 

where if is a detected hadron, typically a pion, [20], and in SIDIS reactions with a trans- 
versely polarized target [21, 22] 

The problem is that in these semi-inclusive reactions ATq a (x) always occurs multiplied by 
a function which depends on the largely unknown Collins fragmentation function. But 
progress is being made and we expect to have a reasonable estimate of A-rg a (x) soon [23]. 

8. SUMMARY 

The standard derivation of the tensorial structure of the expectation value of the angular 
momentum J, for a relativistic spin-s particle, in which the matrix elements of the angular 
momentum operators are related to the matrix elements of the energy-momentum tensor, 
is rendered difficult by the singular nature of the operators involved. The fact that the 
operators are space integrals of densities that contain explicit factors of x^ means that the 
densities do not transform like local operators (we have called these compound operators) 
and has the consequence that the matrix elements of J are highly singular, containing 
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derivatives of delta-functions. Consequently the evaluation of the expectation values of J 
requires a careful limiting procedure, beginning with the off-diagonal elements (p', <?\J\p, cr), 
where a labels the kind of spin state under consideration. We have shown that the results 
in the literature are incorrect, and have derived the correct expressions in three different 
ways, two of them based on a careful wave-packet treatment of the standard approach, and 
the third, quite independent, based on the known rotational properties of the spin states, 
which circumvents the use of the energy-momentum tensor. All three methods yield the 
same results, given in Eqs. (1.1), (1.4) and (1.5). We have emphasized that these matrix 
elements are not Lorentz tensors, but nevertheless the resulting forms have the correct 
Lorentz transformation properties when taken in conjunction with the corresponding forms 
for the matrix elements of the boost operators. We have given some attention to the relation 
between the matrix elements for canonical spin states and helicity states, which is not 
obvious, and requires some care because of the derivatives of the delta functions which 
enter. Using a Fock-space picture of the proton, we have used our results to obtain a new 
sum rule for a transversely polarized nucleon, Eq. (1.10), which involves the transverse spin 
or transversity distribution Arq(x) = h\(x), and which is similar in form to the classic 
longitudinal spin sum rule. 
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